Renormalization Group Study of Chern-Simons Field Coupled to 
Scalar Matter in a Modified BPHZ Subtraction Scheme* 



L. C. de Albuquerque, M. Gomes and A. J. da Silva^ 
Instituto de Fisica, USP 
CP. 66318 - 05389-970, Sao Paulo - SP, Brazil 
(February 1, 2008) 



Abstract 



We apply a soft version of the BPHZ subtraction scheme to the compu- 
tation of two-loop corrections from an Abelian Chern-Simons field coupled to 
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two-loop renormalization group functions are calculated. We compare our 
results with those in the literature. 
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I. INTRODUCTION 



Field theories with the Chern-Simons (CS) term in 2 + 1 dimensions [|IJ,^|] are among the 
best studied models in the past two decades. This is due not only to its potential applications 
but also because of some subtle conceptual and technical aspects. Indeed, the quantization 
of these theories raises interesting questions some of which have not been answered satis- 
factorily up to now. Of particular interest is the set up of a renormalization/regularization 
scheme simple and reliable to provide a consistent framework for complex calculations. Many 
proposals have appeared in the literature 0-0] each of them presenting advantages as well 
disadvantages. For instance, Pauli-Villars regularization explicitly breaks parity, analytic 
regularization is not gauge invariant, and Slavnov regularization becomes rather intricate 
beyond one-loop calculations. 

It is of course desirable that the regularization scheme preserves as much of the models 
symmetries as possible. In this respect, the popular dimensional regularization appears to be 
suited for the task. However, the topological nature of the CS term introduces extra compli- 
cations. The Levi-Civita symbol does not admit a simple extension to complex dimensions 
D and to overcome this problem some modifications have to be done ||. 

If one insists in keeping the Levi-Civita symbol in three dimensions, other approaches 
are possible. Firstly, there is the so called consistent dimensional regularization [[5],|IU| in 
which e^x is treated as essentially three dimensional, but one has to introduce a Maxwell 
or Yang-Mills kinetic term as a supplementary regularization. In another proposal, called 
dimensional reduction fl, [IT|-|i"3f the tensor algebra is done in three dimensions and after- 
wards the Feynman integrals are promoted to D dimensions. This method may introduce 
ambiguities in the finite parts of the amplitudes and also in the divergent parts in high order 
corrections. 

Now, it is possible to renormalize divergent integrals without introducing an specific 
regularization. One of the most efficient and rigorous method is the BPHZ subtraction 
scheme which has been applied in a variety of situations |14]]. One class of this method, the 



soft BPHZ schemes, is specially adequate for the study of massless theories . 

In a soft BPHZ scheme one introduces subtraction operators in the external momenta and 
in the masses of the theory. In order to have a better control of the infrared divergences, some 
of the subtractions have zero mass whereas others have mass equal to [i (the renormalization 
scale). In D = 4 this scheme is in many circumstances equivalent to the dimensional 
regularization with a minimal subtraction prescription, and leads to a mass-independent 
renormalization group equation flIEfl . 

In this work we use a soft BPHZ scheme in an Abelian CS theory coupled with scalar 
matter to compute two-loop renormalization group functions. Since this method do not 
involve analytic continuation in the space-time dimension (i.e., we stay in the physical di- 
mension D = 3), we evade the problems aforementioned. Although we will be dealing mainly 
with massive particles the scheme allows, for non-exceptional momenta, a smooth zero mass 
limit if all super-renormalizable interactions are deleted. The paper is organized as follows: 
In section II we introduce the model, generalize the soft BPHZ scheme to the CS theory 
with scalar matter, and show that the soft BPHZ scheme respects the Ward identity. In 
Section III we compute the renormalization group functions up to two-loops. Finally, we 
draw some conclusions and comments. Our results for the renormalization group functions 
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extend the ones recently computed through the use of the dimensional reduction [13|, in the 



sense that the results agree whenever a comparison is possible. The mentioned agreement 
indicates that dimensional reduction is a consistent scheme, at least up two-loops, as far 
as renormalization group functions are concerned. However, differently from [|L| the model 
that we study includes the most general, renormalizable and U(l) invariant, self-interaction 
of the scalar particles. 



II. THE MODEL AND THE SOFT BPHZ SCHEME 



^012 



The Lagrangian describing the model reads (the metric has signature (+, 
1) 



and 



C = (D^D^ - m 2 ^ - -($ t $) 2 



(1) 



where _D M = d p — ieA^ is the covariant derivative and F pu = A v — d v A^. The canonical 
operator dimensions are as follows (in unity of mass: [m] = 1): [$] = 1/2, [A^] = [A] = [£] = 
1, [u] — [e] — 0. The ultraviolet degree of superficial divergence of a graph T is 



d(T) = 3 - - -N+ - N A , 



(2) 



where is the number of insertions of the super-renormalizable vertex A($^$) 2 , and 
{Na) is the number of scalar (vector) external lines in T. Nonetheless, we must stress that, 
due to the contractions of indexes, this degree is lowered by the number of external trilinear 
vertices having attached one internal CS line. In this work we will use the Landau gauge, 
formally obtained by letting £ — > 0. 

In a soft BPHZ approach the divergent quantum amplitudes, described by Feynman 
integrals in perturbation theory, are made finite by the application of subtraction operators 
T d{~f) arrari g ec [ j n accord with the forest formula 



The subtraction operator r n has degree n in the derivatives and uses derivatives with 
respect to m 2 (each of them having degree 2) if the integrand depends only on the square of 
this mass. All infrared finite subtractions are made at m 2 = but in the would be infrared 
divergent terms (possibly occurring in the last subtractions) we make m 2 = /i 2 . As a rule, we 
make the minimum number of subtractions necessary to render finite a Feynman amplitude. 
Explicitly, the first three subtraction operators are defined as 



T°g(pi,m) = 9(0, fi), 



eg 



r 2 g( Pt ,m) = g(0,0)+p^ 

dg 



T 3 g{ Pi ,m) = g(o,o)+ P i 



Pi=0,m=fj, 

i „ v o 2 g 

Pz=0 ,m=o + 2 PlPj drfdtf 
, 1 , » Q2 G 



d p f Pl=0 ,m=o 2 rtr3 dp f *dp 



(3) 



M 2dg 

+m 7T^ 

Pi =0,m=fi (jm z 

M 2 dg 

+m - — - 

Pi=0,m=0 om z 



Pi=0,m=/i 



Pi=0,m=0 
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i „ „ x d*g 

6 PiPjPk dp?dp'jdp x k 



M 2 „ d 2 g 

+m Pi a Ma 2 
Pi=o,m=ii op^om 1 



Pi=0,m=fi 



It is apparent that our scheme does not introduce infrared divergences in graphs having at 
least one internal scalar line. For ultraviolet divergent graphs without internal scalar lines, a 
special care has to be taken. To be precise, for ultraviolet divergent graphs without internal 
scalar lines, we will use the regularized CS field propagator 



\u(p) 



p 2 - M 2 + ie 



(4) 



where M is a regularization mass which, in the subtraction terms, is to be treated in the 
same way as m and has to be put equal to zero after all the subtractions are performed. 

To see how the method works, in the following we will compute the vacuum polarization 
at one- loop in the theory specified by Eq. (jH). The vacuum polarization is given by the sum 
of two terms, namely: 

n^ = 2ie 2 g pX J[dq]A(q), (5) 
Kfffr) = -e 2 / [dq] (2q + p) p (2q + p) x A(q) A(q + p) = -e 2 J [dq] I pX (p, q; m), (6) 



where A(k) = i/(k 2 — m 2 + ie) is the scalar propagator, and [dq] = d 3 g/(27r) 3 . Both tx 



( B) 

and Ti pX are linearly divergents. For the renormalized quantities we obtain 
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(A) 
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(B) 
pX 



n 



Jrl 
(P) 



j[dq] 


i r 






q 2 — m 2 q 2 



e 2 J [dq](l - r 1 ) I pX {p, q\ m) 
(2q + p) p (2q + p) x 



[dq] 



-4 



QpQx 



(q 2 — m 2 )((q + p) 2 — m 2 ) (q 2 ) 



P x 



dp c 



m=p,p=0- 



(7) 



(8) 



From a practical standpoint, the finite integrals in Eqs. (0) and ® may be most easily 
computed using an intermediate regularization (cutoff, dimensional, Pauli-Villars, etc), the 
final result being of course independent of the regularization employed. Note that the last 
term of the Eq. (@) vanishes upon symmetric integration. We obtain, 



7T 



pX 



e 2 m 



R 



R 



2tx 
8tt 



9pX, 

o ( , VpVx\ P 2 - 4m : 
2m\g pX A ■ 



p- 



sinh 1 ( -j=^ 

\ , /Arn2 



- \/ Am 2 — p r " 



pX 



(9) 
(10) 



where T pX = g pX - p p p x /p 2 . Using that n pX (p) = tt^ + 7i pX '{p), we finally get (note that 
ti {p 2 — >• 0) = as it should) 



IB), 



K pX (p) 



R 



<P)T PX = -i- 



p 2 — Am 2 
2m H ^= — smh 



v 7 4m 2 — p 5 



X, 



P x- 



(11) 
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The Lagrangian density in Eq. ([I]) is invariant under global U(l) transformations. At the 

classical level, this demands the conservation of the current J^x) = £e[<&' $] +2e 2 $^$A /J . 
At the quantum level this current will be quantized with normal product of minimum degree, 
i. e. two. The corresponding Ward identity for the Green functions reads 



d,(TN 2 [J^](x)X {U]} (x k , yj )) = (TN^J^X^x^y,)) 

- N 2N 

J2s{x-Xj) - s ( x 

.j=l j=N+l 



(TX {u . } (x k ,yj)) 



(12) 



where X {u .y(x k ,yj) = Uk=i ®(xk) Ili=N+i&( x i) Yli=iA Ul (yi), follows from the use of the 
normal product algorithm which turns out to be valid in our scheme. In the first step 



Lowenstein's differentiation rule |T^] was used to get the partial derivative inside the normal 
product (which can be done if one increases the degree of the normal product by one unity); 
in the second step one uses the equations of motion for the bilinears <3>tn$ and which 
in our case have the classical form. 

The Ward identity for the proper vertex functions has the same form as in Eq. flT2"|) 
except for a minus sign on its right hand side. Denoting by T^ 2N ' L \pi . . . ,P2N] fci, • • ■ &l) the 
proper vertex function of equal numbers (JV) of $ and <&t fields and L gauge fields, one can 
verify that, in momentum space 

q^^(p,p';q) = -e[T^(p') " T (2) (p)], (13) 

r£ 2 W; i, k) = -e\r?»(p + q, P '; k) - r^(p,p' - 9; k)\. (14) 

Here and in the following we adopt the simplified notation = r^' -*. 



III. RENORMALIZATION GROUP FUNCTIONS 



The renormalized vertex functions introduced in the previous section satisfy the renor- 
malization group equation 



9 a 



d n d n d 1 _ d 



T (N,L) = Q 



(15) 



where @ m 2, (3\, f3 u , /3 e 2, 7$ and ja are power series in the coupling constants e 2 , u, and A. 
We note that j3 e 



as a consequence of the Coleman-Hill theorem [17 



We shall fix now the other functions. We begin by proving that they do not have one-loop 
contributions. Indeed, in the computation of the renormalization group functions the rele- 
vant contributions come through the n dependence of the subtraction terms. As mentioned 
before, these subtractions are those which are potentially infrared logarithmically divergent. 
We shall now examine these possible contributions. We will use a graphical notation in 
which each diagram represents a set of Feynman graphs differing by the orientation of the 
external lines. Moreover, to facilitate the discussion, one may use an auxiliary regularization 
so that each subtraction can be analyzed individually. 

The divergent graphs contributing to I^ 2 ) at one-loop are shown in Fig. ^ Diagram 
|l]a is linearly divergent but, as it does not depend on the external momentum, it does not 
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lead to a logarithmic term. The graphs in |I]b and [T|c vanish upon contraction of indexes. 
More generally, the subtractions for any one-loop graph contributing to T^, with an odd 
number of CS lines, vanishes in the Landau gauge. This can be used to eliminate the possible 
contributions to r( 6 ) coming from Fig. |2|. Finally, Fig. ^ shows divergent contributions to 
at one-loop. Graph |3]a is linearly divergent, with no logarithmic corrections. Graph |3]b 
is linearly divergent. Using in Eq. (|J), and, in accord with @), taking the limit M —>■ 
at the end, we found that it does not depend on Diagram |3]c is actually finite due to the 
observations made after Eq. (0). A similar analysis can be done for r^ 2,1 ) and r^ 2 ' 2 ) leading 
to the result mentioned in the beginning of this section. 

Let be the self-energy function defined by T^(p) = i[p 2 — m 2 — S(p)]. We have 



9 a 



d 



dm 2 



d_ 



+ Px^r +Pu^r -27* 



d_ 

dv 



S(p)+/3 m2 +2 74> (p 2 -m 2 ) = 0. 



It is then easily verified that 



B [1] 



+ 2 7 i 1] (p 2 -m 2 ) 



7$ 



[i] 



4 p 2 
2m 2 7 i 1] 



2 



dp^dp 

d 




p=0 



, yW 



(16) 

(17) 
(18) 

(19) 



p=0 



where the superscript in parenthesis designates the order of the corresponding quantity in 

0, as there is no 



[i] 



[1] 



the loop expansion. From Eqs. ( |T8l) and ([H]) we obtain that 7$ 
H dependence at one-loop. From Eq. (|1^) follows then that ftfi = ffi = 0. 

Let us now proceed to the two-loop calculation of the renormalization group functions. It 
can be easily verified that 7$ and fij 2 are given by Eqs. (|Iq) and (|19D , with the superscript 
[1] replaced by [2]. Besides that, using the one-loop results, and writing the renormalization 
group equations for r( 6 ) and T^ 4 ) in a loop expansion, we obtain at two-loops the following 
relations: 

3|. 2| = 6^| 1 +v|-r <6,[21 (0) (20) 



/3i 21 



6z/7$ ] + , 

4A 7 l 21 +^#r( 4 )[ 2 ](o). 



(21) 



We have used the analytical continuation in the number of dimensions as intermediate 
regularization. It turns out that all integrals needed for the computations may be expressed 
in terms of (n-dimensional Euclidean space, with {dq} = a € d n q/(2ir) n , where a is a mass 
scale, and e = 3 — n) 



l[ n) (M a ,M b ,M c 



{dk}{dq} 



1 



(k 2 + M 2 )(q 2 + M 2 )((k + q) 2 + M 2 



32tt 2 



- -7 + 1 
e 



In I 



(M a + M b + M ( 
Ana 2 




J 2 (n) (M a ,M 6 )= f{dk}{dq} 



(k 2 + M 2 )((k + q) 



Ml 



I67P 



M a M b . 



(22) 
(23) 
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At two-loops, the nonvanishing contributions to 7$ come from the three (quadratically 
divergent) diagrams in Fig. [|. The diagram |4]a may be written as 



£g ] (p) = -2te 4 J [dk] [dq] A(k + q- p) A^(k) A^(q) 



(24) 



SBPHZ 



This diagram has three divergent subgraphs and eight forests, but the relevant term 
comes from the subtraction associated to the diagram as a whole. An application of the 

rnl 

forest formula leads to the following contribution to 7$ 

„i 



[2] 
7$ A 



48tt 2 



(25) 



The second graph, Fig. £|b, corresponds to the unsubtracted integral 

S| ] (p)= ie A J [dk}[dq] (2p + ky{2p + 2k + q) a (2p + k + 2q) v {2p + qf A(p + k) 



xA(p + k + q)A(p + q) A^{k) A af3 {q) 



(26) 



SBPHZ 



In this case there is no divergent subgraph and the calculation of 7$ B is reduced to just the 

(27) 



contribution from the forest corresponding to the graph as a whole. We then get 

„i 



[2] 

70s 



12tt 2 



The graph in Fig. ^c is written as 

Eg 1 (p)= ie 4 J [dk] [dq] (2p - ky(2q - k) a {2q - kf{2p - k) u A(k - p) 
xA(q)A{q-k)A^ a (k)A f3u {k) 



(28) 



SBPHZ 



The calculation now is more involved and details will provided in the Appendix. That 
analysis produces the result 



[2] 

7*c 



24vr 



2 ' 



Adding the results in Eqs. 



[2] 

7$ 



and (p9|) , we obtain 
7e 4 



48tt 2 



(29) 



(30) 



which is in accord with the result obtained in |13[ using dimensional reduction. 

The diagrams that contribute with a logarithmic correction to r^ 6 ^ 2 ' are shown in Fig. 
[5]; each one of them is superficially logarithmically divergent. Concerning these graphs we 
make the following comments. Diagrams |5]a and |5]d do not have divergent subgraphs. The 
other graphs, [5|b-c and [|e-f have just one divergent subgraph. However, in all cases the 
contribution coming from the subgraph vanishes. For the graph ||b this happens after the 
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contraction of indexes whereas for the graphs [|e and [|f this results from a cancellation 
between different forests; also, the possible contribution arising from the subgraph of |5|c 
is just the graph |3|b whose // dependent subtractions vanish under symmetric integration. 
The conclusion is that in each case one has to compute only the contribution of the forest 
containing just the graph as a whole. The results of these calculations are summarized in 
table |. Summing those contributions and using Eq. (|2Cf) we obtain 



This result coincides with that of [13]. To obtain (3^ we need to compute logarithmic 
contributions arising from r^^. The relevant graphs are listed in Fig. |[ The graph |6|a has 
d{T) = 0, and does not contain a divergent subdiagram. A straightforward calculation gives 



^-^"H^J- (32) 

The diagram []b also has d(T) = 0, and no divergent subdiagram. It leads to 

r<*' = i^lng). (33) 

The diagram ||c has the diagram |3]b as a divergent subdiagram. We already know that this 
subgraph does not contribute. The calculation of the the forest containing just the overall 
diagram gives 

r ™ = _i_ Ae Mn(^). (34) 
Taking the above results and using Eq. (^1|) we obtain 

We turn now to the computation of (3^j 2 - According to Eq. fliTf), we still need to compute 
the graphs shown in Fig. [7|. The diagram [7|a is the same as that in Fig. |]a, now with p = 0. 
There is a cancellation between forests. After some manipulations, we obtain 

£2 ] (0)= * m 2 eMn(^). (36) 

The diagram in 0b has d{T) = 2, and three divergent subgraphs. However, due to cancella- 
tions the end result is 

Ef(0) = - i i 5 mV !„(=). (37) 
Finally, graph 0c has d(T) = and no divergent subdiagram. The result is 

^°) = 3^ 2 0- (38) 



S 



Collecting the partial results, £^(0) = sf(0) + s| ] (0) + sg ] (0) from Eqs. 
and using Eq. ([0]) with n = 2, we obtain 



For A = this agrees with the anomalous dimension of the composite operator $t<I> as 



computed in |]13| , as it should. 

To discuss the fixed points structure of the model we introduce a dimensionless coupling 
by A = A/i. Up to two- loops the beta functions are given by Eq. (|3ll) , 



A = -A + ^-^e*A, (40) 

and by Eq. (^) with A replaced by /iA. If m and A are zero there are no induction of 
and ($^<3>) 2 counterterms as (3^ and f3 m 2 both vanish. The case e = has been analyzed in 
the literature |18|] unveiling an interesting tricritical behavior. Near the trivial fixed point, 
for small momenta the running couplings A e / and v e f are driven away and approach the 
origin, respectively. For large momenta the opposite happens. As mentioned in ||13|| , for 
e ^ there are no other fixed points since, f3 v never vanishes in the perturbative region. 
Besides, the behavior of X e f near the origin is not sensible to the introduction of e. 



IV. CONCLUSIONS 

In this paper we have shown that it is possible to define a consistent, gauge invariant 
subtraction scheme with a soft behavior in the infrared regime (soft BPHZ) for an Abelian 
Chern-Simons theory coupled to a massive scalar matter with A($t$) 2 and z/($"l"$) 3 self- 
interactions. Within the soft BPHZ approach we circumvent the problems associated with 
the analytical continuation of the Levi-Civita tensor. Hence, there is no need to deal with 
a complicated consistent dimensional regularization and neither it is necessary to introduce 
a Maxwell term. However, there is a price for this simplification, since we have to deal 
with various forests with the consequent increase in the number of Feynman integrals. In 
the process of calculation, however, we have found that these integrals can be reduced to a 
few primitive ones. In all cases that we studied we explicitly verified the finiteness of the 
subtracted integrals. 

We have done a two-loop calculation of the renormalization group functions. Analogous 
models were studied in |||TJ,[I3|,[n|. Our 7$' agrees with the one computed in || (Abelian 



case). We found only a qualitative agreement with |TT[ but a more close comparison seems 



infeasible due to the lack of details in 11 



The comparison with the renormalization group functions computed with the consistent 



dimensional regularization scheme in [|19[] is more difficult. Indeed, the RG functions com- 
puted in |19[ contain divergent contributions in the pure CS limit (no Maxwell term) that 
are absent in our approach. However there is some partial agreement between our results 
and the finite parts of /3^\ f3^ 2 , and /3ffl of [|19|] (some coefficients of the expansion of these 

[21 



functions are identical to ours). Our result for 7^ is entirely different from that in [|19 



as 



2 

those authors claim that 7$ = 0. The discrepancy could in principle be attributed to the 
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use of different renormalization schemes. The use of an extra regularization represented by 
a Maxwell term for the field brings additional complications in their proposal as some of 
the coefficients of the renormalization group functions become singular as the regularization 
is removed. A more careful analysis of this method is still lacking. 

Although our main interest resides in the pure CS model, we would like to make a few 
comments on the possible changes in our scheme if a Maxwell term —jF^F^ is added to 
Eq. (0). The propagator for the gauge field is then modified to 

pp i 

V(p) = W p2 _ M 2 + l- a 2p2 ( a2 W>" + iaT ^)> ( 41 ) 

where M is the auxiliary massive parameter as introduced in Eq. (^) and T^ v is the transver- 
sal projector defined below Eq. (|T0|) . Observe that a has dimension —1 in units of mass. 
The degree of superficial divergence is now given by 

d{Y) =3-^N A - iiV - - V AH , - l -V A ^ (42) 

where Vq denotes the number of vertices associated to O in T. With this new power counting 
many graphs previously divergent turn out to be convergent. In this situation, the following 
possibilities can be envisaged. 

i) Instead of (0) one adopts (|4~2"| ) to define the graphs to be subtracted. The outcome is 
a well defined theory as far as a, the coefficient of the Maxwell term, is kept nonvanishing. 
The result however is not analytic in a and the limit a — > does not exist. 

ii) One uses the old power counting and the parameter a is not changed in the subtraction 
terms. This means that many graphs will be oversubtracted. However such subtractions are 
needed if, as a — ► 0, one wants to recover the pure CS model studied in this paper. 

Finite one-loop renormalization constants for the non-Abelian CS theory with fermionic 
matter were computed in j2(J using consistent dimensional renormalization and found to be 
different from the ones obtained from the dimensional reduction prescription. It should be 
noticed that, as remarked in [J21]J , even for finite theories there could exist different fami- 
lies of BRST invariant regularizations leading to distinct results for the one-loop radiative 
corrections. This may imply that in spite of the numerical differences there are no physical 
inconsistencies between the two approaches. 

The model studied in [[13|] is a particular case of Eq. ([!]), with A = m 2 = 0. We found a 
complete agreement in 7$, (3]?\ and 7$j $ . O ur results show that the dimensional reduction 
method is consistent at two-loops, at least, as far as the computation of the RG functions 
is concerned. 
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APPENDIX A: 



In this appendix we shall present details of the calculation of the contribution of the 
diagram ^c to 7$'. The graph (T) is logarithmically divergent and has just one divergent 
subgraph which is the same as the one associated to the 7r^ contribution to the CS vacuum 
polarization. We will denote this subgraph by 7. There are four forests: (the empty 
forest), T, 7 and {V, 7}. Thus, the application of the forest formula to (p8|) leads to 

sg ] (p) = 16ieVp" / \ dk \ [dq}{q a q p A fia (k)A (Sl/ (k) [A(k - p, m)A(g, m)A(q - k, m) 
—A(k, n)A(q, Li)A(q - k, fx)] - q a q l3 A fia (k)A Pu (k)A(k - p, m) [A 2 (g, 0) - ik ■ qA 3 (q, fx) 
+q a q^A fUX {k) A pv {k) A{k, fx)) [A 2 (q, 0) - ik ■ qA 3 (q, fx)] } (Al) 

The divergent parts of the terms containing the k ■ q factor cancel between themselves; the 
finite parts associated to them are odd and vanish under symmetrical integration. We can 
therefore rewrite the above expression as 



Y%\p) = lGeWe^e^ J [dk][dq]k x k^ 



1 



[{k — p) 2 — m?\ [q 2 — m 2 ] [(q — k) 2 — m 2 ] 



[k 2 - fi 2 ][q 2 - fi 2 }[{q - k) 2 - fx 2 ] {{k - p) 2 - m 2 ]{q 2 ) 2 \(k 2 - fx 2 }{q 2 ) 

[21 

The contribution to 7$ arising from the above integral can be most easily computed if 
one employs an intermediate auxiliary regularization. Adopting dimensional regularization 
we arrive at 

7 i| = -16e 4 [X 1 -X 2 ], (A3) 

where 

Xl = /[^][^| I fc2(fcV)(g2 _/ 2)((fc + g)2 _ /i2r (A4) 

(k 2 ) 2 (k 2 - fx 2 )(q 2 - fx 2 )((k + q) 2 - fx 2 ) (A5) 

which, after some simple manipulations, can be calculated with the help of the results 
and d|). We obtain 

XAu) = (A6) 
w 960vr 2 ' v ; 

X 2 (fx) = l — , (A7) 

vw 640tt 2 ' v ; 

from which we obtain the result quoted in the text. 
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TABLES 

TABLE I. Contributions to L^ arising from the graphs in Fig. ||| 

Diagram type Contribution 

5b £ In (f 

5c M e 4 u In ( ^ 

5d -^eMnf^ 

5e -4e 8 ln(^ 

5f -4e 8 ln( 
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FIGURES 




a b c 

FIG. 1. r( 2 ) at one loop. 




a b 
FIG. 2. Divergent contributions to at one loop. 




a b c 

FIG. 3. Divergent contributions to at one loop. 




a b c 

FIG. 4. Contributions to 7$. 
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d e 

FIG. 5. Two-loop contributions to T^ 1 . 




a be 
FIG. 7. Contributions to Y^{p = 0). 
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